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I.  Introduction 


1.  Background 

Helicopters  are  frontline  systems  in  most  of  the  world’s  armed  forces,  with  combat  mis¬ 
sions  including  troop  insertion,  observation  and  attack.  Operational  tactics  typically  employ 
low  level  flight  to  mask  the  helicopter’s  presence.  Although  this  can  conceal  the  aircraft 
from  long  range  threats,  the  aircraft  may  still  be  exposed  to  attack  from  nearby  small  arms 
and  light  anti-aircraft  artillery  (e.g.,  23  and  30  mm).  Combat  experience  shows  that  rotor 
blades,  particularly  those  of  the  main  rotor,  are  often  hit.  While  the  rigorous  application  of 
damage-tolerant  materials  and  vulnerability  reduction  design  principles  yields  high  levels  of 
survivability,  ballistic  damage  to  the  blade(s)  (particularly  from  high-explosive  projectiles) 
may  still  affect  the  performance  capabilities  and  airworthiness  of  the  vehicle. 

Vibrations  generated  by  the  damaged  rotor  system  may  be  severe  enough  to  compro¬ 
mise  the  structural  integrity  of  rotor/hub  subsystem,  to  cause  the  crew  to  abort  the  mission, 
or  force  the  aircraft  to  land.  Damage  to  the  rotor  subsystem  may  also  induce  an  aeroelastic 
instability  that  can  cause  loss  of  flight  control  and  an  attrition  kill.  Currently,  with  resp>ect  to 
vulnerability  assessment,  there  is  insufficient  understanding  of  rotor  blade  ballistic  damage 
effects  on  the  aeroelastic  response/stability  of  the  rotor  system  and  the  associated  helicopter 
performance  and  vibratory  loads.  An  understanding  of  these  effects  is  not  only  essential  to 
accurately  assess  the  level  of  vulnerability  but  also  benefits  helicopter  survivability  improve¬ 
ment  efforts. 

Most  previous  ballistic  vulnerability  studies  involving  helicopter  rotors  were  simply 
conducted  on  isolated  blades,  which  were  statically  evaluated,  sometimes  with  axial  loading. 
Using  thic  approach,  one  can  locally  identify  the  effects  of  damage  on  the  blade  structural 
properties,  such  as  the  changes  in  local  strength  and  strain  or  residual  stress.  Those  studies, 
however,  may  not  be  directly  applied  to  assessing  the  performance  of  a  whole  helicopter 
(i.e.,  with  rotor-fuselage  coupling)  suffering  blade  damage  in  flight.  In  the  present  study, 
a  specialized  aeroelcistic  analysis  has  been  performed  to  investigate  helicopter  response  to 
in-flight  main  rotor  blade  damage. 

In  the  context  of  the  Army  Research  Laboratory’s  (ARL)  Process  Structure  for  analyz¬ 
ing  combat  system  (e.g.,  helicopter)  vulnerability,  this  research,  and  its  associated  engineer¬ 
ing  based  methods,  ciddresses  the  mapping  from  Level  2,  the  Target  Component  Damage 
State  to  Level  3,  the  Target  Capability  State  (i.e,  mapping);  other  levels  include  1, 
Threat/Target  Initial  Conditions  and  4,  Target  Combat  Utility.  A  significant  feature  of 
this  process  is  that  at  each  level  (or  space)  distinct,  measurable  information  is  available 
defining  the  threat/target  encounter  and  vulnerability /lethality  outcome.  For  example,  here 
physical  and  aerodynamic  factors  defining  rotor  blade  damage  (Level  2)  are  mapped  via  en¬ 
gineering  methods  into  parameters  which  define  the  rotor  and  helicopter  system’s  functional 
capability  (Level  3);  all  the  defining  terms  are  explicit  and  measurable  through  experimen- 
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tation.  Note  also  that  ambient  (e.g,  air  density)  and  operating  factors  (e.g.,  forward  flight 
velocity)  (defined  at  Level  1)  are  accounted  for  in  the  methodology  and  provide  realistic 
variable  scenario  capability.  Application  of  these  and  other  engineering  analysis  tools  to  the 
Vulnerability/Lethcdity  Process  Structure  will  occur  largely  through  implementation  in  the 
Degraded  States  Vulnerability  Methodology  for  level  02,3  mapping  now  under  development 
for  aircraft  targets  at  ARL. 


2.  Present  Study 


The  objective  of  the  current  research  is  to  investigate  the  effects  of  (simulated)  ballistic 
damage  on  helicopter  rotor  blade  response  and  hub  loads. 


This  investigation  was  performed  using  the  UMARC  helicopter  aeroelcistic  analysis  com¬ 
puter  code  [1].  A  finite  element  procedure  is  used  for  a  structural  analysis  of  damaged  rotor 
blades,  and  rotor  aerodynamic  loads  are  calculated  using  quasisteady  aerodynamic  theory. 
Dynamic  responses  of  helicopter  rotor  subsystems  are  calculated  from  periodic  normal  mode 
equations  using  a  finite  element  method  in  time.  Also,  in  the  present  study,  results  of  wind 
tunnel  tests  of  the  damaged  rotor  blades  [2]  are  used  to  include  the  aerodynamic  performance 
degradation  of  a  ballistically  damaged  blade. 

Results  are  first  calculated  for  the  SA349/2  Gazelle,  an  advanced  geometry  three-bladed 
articulated  rotor  helicopter,  with  undamaged  blade  configuration.  Results  then  are  calcu¬ 
lated  for  this  helicopter  with  simulated  blade  damage  configurations.  The  effects  of  this 
damage  on  helicopter  response  and  loads  are  determined  in  terms  of  blade  flap-lag-torsion 
frequency  and  modal  shapes,  blade  aeroelastic  response,  and  hub  load  variations. 


The  present  study  includes  sever  J  new  features  that  have  not  been  previously  addressed 
in  the  literature. 

1.  Consistent  helicopter  rotor  blade  damage  modeling,  including  eierodynamic  and  structural 
interactions,  is  incorporated  in  the  aeroelastic  analysis. 

2.  Aerodynamic  properties  of  damaged  blade  estimated  from  the  results  of  wind  tunnel  tests 
with  aetual  helicopter  hardware. 

3.  The  ultimate  issues  of  aircraft  flyability  and  controllability  for  given  damages  are  discussed 
in  the  context  of  aircraft  vulnerability  assessment. 
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II.  Formulation 


The  baseline  helicopter  aeroelastic  code  used  in  this  investigation  is  a  comprehensive 
aeroelastic  analysis  based  on  finite  element  theory  in  space  and  time  [1]. 

The  rotor  blade  is  assumed  to  be  an  elastic  beam  undergoing  flap  bending,  lag  bending, 
elastic  twist,  and  axial  deflections.  The  helicopter  fuselage  is  assumed  to  be  a  rigid  body 
undergoing  six  degrees-of- freedom  movement.  Each  blade  is  discretized  into  a  number  of 
beam  elements,  and  for  each  element,  there  is  a  continuity  of  displacement  and  slop>e  for 
flap  {w)  and  lag  (v)  deflections  and  a  continuity  of  displacement  for  axial  (u)  and  torsion 
(^)  deflections.  There  are  two  internal  nodes  for  axial  displacement  and  one  for  elastic 
twist  resulting  in  a  total  of  15  degrees  of  freedom  for  each  element  (see  Figure  1).  The 
formulation  for  the  blade  and  fuselage  equations  of  motion  is  based  on  Hamilton’s  principle. 
The  analysis  is  developed  for  helicopters  with  nonuniform  rotor  blades  having  pretwist, 
precone,  and  chord-wise  offsets  of  the  center  of  mass,  aerodynamic  center  and  tension  center 
from  the  elastic  axis. 


Aerodynamic  loads  are  calculated  using  quasisteady  strip  theory.  Noncirculatory  aero¬ 
dynamic  forces  are  also  included.  To  include  the  effect  of  high  angle  of  attack  flows,  a 
dynamic  stall  model  proposed  by  Johnson  [3]  is  incorporated.  Dynamic  stall  characterizes 
the  delay  in  flow  separation  because  of  unsteawly  angle  of  attack,  and  the  shedding  of  a  vortex 
from  the  leading  edge  of  the  airfoil  when  it  gets  into  a  deep  stall  condition.  These  effects  are 
introduced  in  the  calculation  of  section  lift,  drag  and  pitching  moment.  For  this,  the  time 
history  of  blade  motion  of  the  previous  cycle  is  used. 

Aerodynamic  coefficients  are  computed  in  the  form  of  analytical  expressions  as  well  as 
data  tables.  These  are  represented  eis 


Cl  =  Co(a,A/,q,q)  +  Ci(Q,M,q,q)Q 
Cd  =  </o(a,A/,q,q)  +  ‘/i(Q',A^,q,q)a 

+d2{Q,M,q,q)a‘^  (1) 

Cmac  =  /o(a,A/,q,q)-f /,(Q,M,q,q)Q 


in  which  q  and  q  are  arrays  of  nodal  displacements  and  velocity  vectors,  respectively.  Equa¬ 
tion  (1)  represents  a  set  of  nonlinear  aerodynamic  coefficients  including  effects  of  blade 
motion  and  Mach  number. 
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1.  Coupled  Trim  Analysis 


The  coupled  trim  analysis  consists  of  two  phases,  vehicle  trim  and  steady  response, 
calculated  as  one  coupled  solution  using  a  modified  Newton  method.  For  a  given  flight  con¬ 
dition,  the  control  settings  and  the  blade  steady  response  must  satisfy  both  the  blade  and  the 
vehicle  equilibrium  conditions.  The  method  of  solving  simultaneously  the  blade  responses 
and  the  trim  control  settings  is  referred  to  as  the  coupled  trim  analysis.  An  uncoupled  trim 
solution  based  on  the  rigid  flapping  blade  assumption  is  used  as  an  initial  estimate  for  the 
coupled  trim  analysis.  With  the  trim  control  settings,  the  blade  steady  responses  are  calcu¬ 
lated.  Using  the  blade  responses,  the  rotor  hub  loads  and  a  new  vehicle  equilibrium  position 
are  recomputed.  The  control  settings  are  then  updated  btised  on  the  new  equilibrium  con¬ 
dition. 


2.  Vehicle  Trim 

Propulsive  trim,  which  simulates  an  aircraft  free  flight  condition,  is  used  to  calculate 
the  initial  rotor  control  settings.  The  solution  is  determined  from  the  overall  equilibrium 
equations  (see  Equation  2):  three  force  (vertical  Z,  longitudinal  X  and  lateral  T),  and  three 
moment  (pitch,  roll  and  yaw)  equations.  These  are 


F\  =  Fj^  -f-  Dp  cos  ajfp  —  T  sin  q,  =  0 
Fa  =  F^  +  Yp  cos  <i>,  —  Vir  +  T  sin  =  0 
Fa  =  F^  —  T  cos  o,  cos  —  Dp  cos  sin  qhp 

+  Yp  sin  =  0  (2) 

Fi  =  +  Mxp  -k-T  {h  cos  a,  sin  <f>t  -  ycg  cos 

-h  Yp  cos  a,  cos  -I-  ycg  sin  =  0 
Fj  =  A/^  -I-  Myj.  +  T(Ji  sin  a,  -  Xcg  cos  q,) 

+  Dp  (—ft  cos  Ofip  -h  Xcg  sin  qhp)  =  0 
Fe  =  Mg-\-Mxr-XtYtr  =  0 
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in  which  Fi,F2,  and  F3  are,  respectively,  the  force  equilibrium  equations  in  the  X,Y,  and 
Z  directions,  and  F4,  F5,  and  Fe  are  the  rolling,  pitching,  and  yawing  moment  equilibrium 
equations,  respectively.  Also,  Dp  is  the  fuselage  drag;  Yp  is  the  fuselage  side  force;  Ytr  is  the 
tail  rotor  thrust;  T  is  the  main  rotor  thrust;  Xcg  and  ycg  and  h  are,  respectively,  the  relative 
location  of  the  rotor  hub  center  with  respect  to  the  vehicle  center  of  gravity  in  the  X,  Y,  and 
Z  directions;  Xt  is  the  nondimensional  length  (divided  by  rotor  radius)  between  the  maiin 
rotor  hub  and  the  tail  rotor  hub;  and  a,  and  are  the  longitudinal  and  latertil  shaft  tilts, 
respectively  (see  Figure  2).  Furthermore, 


=  oiffp  -  6pp 


(3) 


in  which  a^p  is  the  hub  plane  tilt  relative  to  the  flight  direction,  and  dpp  is  the  flight  path 
angle  relative  to  the  longitudinal  axis. 


For  a  specified  weight  coefficient  C^,,  the  unknown  quantities  to  be  determined  from 
the  vehicle  equilibrium  equations  are 


VirJ  . 


(4) 


These  values  are  recalculated  iteratively  using  the  modified  hub  forces  and  moments 
including  the  blade  elastic  responses.  The  solution  technique  is  based  on  a  modified  Newton’s 
method.  The  rotor  controls,  which  are  updated  at  the  iteration,  can  be  expressed  as 


Ui+i  =  Ui  +  Aui 


(5) 


in  which 


Aui  =  - 


9u 


F(ui) 


«=«o 


(6) 


in  which  Uo  are  the  trim  control  settings  obtained  initially  using  the  rigid  flapping  blade 
solution,  and  the  Jacobian  5F/5u  is  calculated  using  the  finite  difference  approach.  For 
computational  efficiency,  the  Jacobian  is  computed  only  once  initially  and  used  for  subse¬ 
quent  iterations. 
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3.  Rotor  Dynamic  Response 


The  rotor  dynamic  response  involves  the  determination  of  time-dependent  blade  po¬ 
sitions  at  different  azimuth  locations  for  one  rotor  revolution.  To  reduce  computational 
time,  the  finite  element  equations  are  transformed  into  modal  space  as  a  few  normal  mode 
equations  using  the  coupled  natural  vibration  characteristics  of  the  blade.  These  nonlinear 
periodic  coupled  equations  are  solved  for  steady  response  using  a  finite  element  in  time  pro¬ 
cedure  based  on  Hamilton’s  principle  in  weak  form  [4].  One  rotor  revolution  is  divided  into  a 
number  of  azimuthal  elements,  and  then  periodicity  of  resp>onse  is  used  to  join  the  motions  of 
the  first  and  last  elements.  The  assembly  of  elements  results  in  nonlinear  algebraic  equations 
that  are  solved  using  the  Newton- Raphson  procedure. 

After  the  blade  response  is  obtained,  hub  forces  and  moments  are  calculated  using  the 
force  summation  method.  These  values  are  finally  used  to  recalculate  the  coupled  trim  con¬ 
trol  values. 


4.  Rotor  Wake  Modeling 

For  the  induced  in-flow  distribution  on  the  rotor  disk,  a  free  wake  model  [5]  is  fully 
coupled  in  the  rotor  aeroelastic  analysis.  The  model  can  account  for  wake  self-distortion  by 
updating  its  geometry  according  to  newly  calculated  in-flow  and  blade  circulation  distribu¬ 
tions.  The  geometry  of  the  free  wake  is  divided  into  three  regions:  near  wake,  rolling  up 
wake,  and  far  wake.  The  near  wake  consists  of  a  series  of  radial  panels,  each  with  linear 
circulation  distributions.  The  rolling  up  wake  consists  of  an  in-board  linear  circulation  dis¬ 
tribution  panel,  and  a  tip  panel  that  represents  the  rolling  up  of  the  tip  vortex.  The  far  wake 
is  modeled  as  one  panel  with  a  linear  circulation  distribution,  and  a  concentrated  tip  vortex 
whose  strength  is  proportional  to  the  maximum  circulation  value  on  the  rotor  blade.  The 
helical  geometry  of  the  concentrated  tip  vortex  is  updated  while  the  in-board  wake  portions 
are  not  changed. 

The  free  wake  analysis  is  implemented  in  three  stages.  First,  blade  motion  and  loading 
are  calculated  using  a  linear  in-flow  model.  Next,  wake-induced  coefficients  are  calculated 
for  an  undistorted  wake  geometry.  The  nonuniform  in-flow  is  calculated  and  is  used  to 
obtain  blade  motion  and  loading.  Finally,  the  free  wake  geometry  is  calculated.  For  this,  the 
influence  coefficients  are  re-evaluated,  and  blade  motion  and  loading  are  again  obtained  using 
nonuniform  in-flow  values.  For  subsequent  iterations,  the  free  wake  geometry  is  generally 
held  fixed,  and  only  the  strength  of  vortices  is  updated. 


III.  Solution  Procedure 


The  following  procedure  was  used  to  study  the  simulated  ballistic  damaged  effects  on 
helicopter  rotor  response  and  loads  in  forward  flight  conditions. 
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1.  With  prescribed  input  data,  vehicle  trim  equations  are  calculated  using  rigid  blade 
flapping  (a  starting  point). 

2.  Using  control  inputs  from  the  vehicle  trim  solution  of  Step  1  and  prescribed  blade 
damage  condition,  the  blade  nonlinear  steady  response  is  calculated.  The  results  give  de¬ 
tailed  individual  blade  responses  at  different  span-wise  and  azimuthal  positions. 

3.  Hub  loads  and  moments  are  calculated  using  elastic  rotor  responses.  Then,  the 
vehicle  trim  values  and  blade  responses  are  recalculated  iteratively  using  the  modified  hub 
forces  and  moments.  This  step  is  repeated  until  a  converged  solution  is  obtained. 

4.  Steps  of  2  to  3  are  repeated  for  different  damage  configurations. 


IV.  Results  and  Discussion 


Numerical  results  are  first  calculated  for  a  undamaged  (baseline)  rotor  of  the  SA349/2 
Gazelle  helicopter  for  level  flight  conditions.  Results  then  are  calculated  for  this  rotor  with 
simulated  damaged  conditions  and  the  effects  of  blade  damages  are  assessed.  Some  of  im¬ 
portant  structural  and  aerodynamic  characteristics  of  this  helicopter  rotor  are  given  in  Table 
1.  Other  characteristics  of  this  helicopter  are  given  in  Reference  6. 

The  SA349/2  Gayelle  helicopter  wais  fitted  with  an  advanced  geometry  three-bladed 
articulated  rotor.  The  Grande  Vitesse  (GV)  blades  consisted  of  OA209  advanced  airfoils 
(9%  thickness  ratio),  adjustable  tips,  and  nonlinear  twist  distributions. 

In  the  present  study,  an  effort  wcis  meide  to  model  the  helicopter  blade  characteristics 
as  accurately  as  possible.  The  structural  properties  of  blades  obtained  from  the  measured 
data  [6]  were  used  for  blade  finite  element  model  (see  Table  2).  Nonlinear  twist  distribution 
of  the  blade  was  also  modeled  in  the  blade  finite  elements.  One  of  the  main  features  in  this 
articulated  rotor  blade  is  a  lag  damper.  The  lag  damper  can  significantly  affect  the  blade 
dynamics.  In  the  present  analysis,  additional  stiffness  as  well  as  damping  t^rms  attributable 
to  the  lag  damper  were  included. 

For  the  calculation  of  blade  dynamic  response,  each  rotor  blade  is  discretized  into  eight 
beam  elements,  and  each  beam  element  consists  of  15  nodal  degrees  of  freedom.  For  normal 
mode  reduction,  five  coupled  rotating  natural  modes  are  used,  comprised  of  tw’o  flaps,  two 
lags,  and  one  torsion  mode.  For  periodic  response,  one  cycle  of  time  is  discretized  into  eight 
time  elements,  and  each  time  element  represents  a  fifth  order  ix)lynomial  distribution  of 
motion. 


Table  1:  SA349/2  Gaselle  Helicopter  Characteristics. 


Aircraft  gross  weight 

4400  lbs 

Number  of  blades,  Nf, 

3 

Radius,  R 

17.5  ft 

Blade  chord,  c 

1.17  ft 

Solidity,  a 

0.064 

Lock  number,  7 

5.13 

Rotor  thrust  ratio,  Crjcr 

0.064 

Blade  airfoil 

OA209 

Rotational  speed,  Q 

387  rpm 

Rotating  flap  natural  frequency 

1 .02/rev 

Rotating  lag  natural  frequency 

0.49/rev 

Rotating  torsion  natural  frequency 

4.28/rev 

Fuselage  lift  coefficient,  Ctf 

0.00499 

Fuselage  drag  coefficient,  Cof 

0.0092 

Fuselage  pitch  moment ,  Cm/ 

-0.00086 

Fuselage  roll  moment ,  Cij 

-0.0001 

Fuselage  yaw  moment ,  Cn, 

0.0001 

Table  2:  Structural  Properties  of  Gazelle  Rotor  Blade. 


1.  Baseline  (Undamaged)  Helicopter:  Correlation  Study 


To  validate  the  coupled  trim  procedure  used  in  the  present  study,  calculated  trim  values 
of  baseline  (undamaged)  helicopter  are  correlated  with  existing  flight  test  data  [6]  for  four 
level  flight  conditions  (advance  ratio  ^  =  0.14,  0.256,  0.333,  0.378). 


Figure  3  shows  the  correlation  of  the  power  requirement  to  maintain  the  level  flight  of 
the  helicopter.  The  torque  coefficient  {Cq)  obtained  from  the  coupled  trim  solution  was  used 
to  calculate  shaft  power.  The  quasi-steady  aerodynamics  with  the  free  wake  model  are  used 
for  the  entire  flight  envelope.  Excellent  correlations  with  flight  test  data  are  observed. 


The  blatde  flapping  angles  are  also  correlated  with  flight  data  and  shown  in  Figures  4 
through  6.  Calculated  blade  flap  responses  at  the  tip,  in  terms  of  the  coning  angle  (^) 
and  the  rotor  disk  tilt  angles  {0ic,  0ii),  are  compared  with  the  data  for  four  flight  condi¬ 
tions.  Reasonably  good  agreement  with  test  data  was  observed,  providing  the  author  with 
a  certain  confidence  in  the  helicopter  rotor  subsystem  model  constructed  in  the  present  study. 


In  the  following  section,  results  are  calculated  for  the  helicopter  damaged  on  its  rotor 
blade,  and  then  blade  damage  effects  on  blade  and  hub  loads  are  assessed. 


2.  Helicopter  with  Damaged  Blade 


Blade  ballistic  damage  is  simulated  as  a  reduction  in  structural  properties  as  well  as 
a  change  in  aerodynamic  characteristics  of  blade.  The  structural  damage  is  simulated  as  a 
50%  reduction  in  bending  and  torsional  stiffness  {EIy,EIt,  and  GJ)  and  mass  of  damaged 
blade  section.  The  aerodynamic  damage  is  simulated  as  a  change  in  blade  aerodynamic 
characteristics  in  terms  of  section  lift,  drag,  and  moment  coefficients.  In  the  present  study, 
the  following  relations  are  used: 


G/damage 

^(^damage 

damage 


=  01  G/ undamaged 
—  0d  ^(^undamaged 
0m  Gm^c  undamaged 


(7) 


in  which  0u0di0m  are  respectively  the  lift,  drag,  and  moment  deficiency  functions.  The 
values  used  in  the  present  study  are  summarized  in  Table  3.  These  are  estimated  based  on 
the  wind  tunnel  test  results.  The  details  of  test,  including  data  acquisition  and  reduction, 
arc  available  in  Reference  2. 
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Table  3:  Aerodynamic  Deficiency  Functions  of  Damaged  Airfoil. 


Deficiency  Functions 

A/.nc  <  0.7 

Mine  >  0.7 

A 

0.75 

2.5 

0m 

0.9 

0.9 

TTTmJ 

Table  4:  Frequencies  of  Damaged  Blades  (per  rev). 


Mode 

Baseline 

Damage  1 

Damage  2 

First  Lag 

0.492 

0.575 

0.856 

First  Flap 

1.016 

1.048 

1.181 

Second  Flap 

2.694 

2.783 

2.923 

First  Torsion 

4.289 

4.460 

4.438 

Second  Lag 

4.611 

4.768 

5.056 

In  the  present  study,  two  different  blade  damage  conditions  are  considered;  Damage 
1  and  Damage  2.  Damage  1  represents  a  condition  in  which  10%  of  the  outboard  blade  tip 
of  one  blade  is  damaged,  resulting  the  reduction  in  structural  and  aerodynamic  capability 
of  that  blaxle  section.  Damage  2  represents  a  more  severe  damage  condition,  in  which  20% 
of  the  outboard  section  is  damaged  (see  Figure  7). 


Effects  of  the  simulated  damages  on  the  bI^uie  natural  frequencies  are  summarized  in 
Table  4,  and  the  corresponding  mode  shapes  are  shown  in  Figures  8  through  12.  The  results 
of  baseline  (undamaged)  blade  are  shown,  along  with  two  damage  cases.  It  is  shown  that  the 
blade  damage  affects  all  the  modes  of  the  blades.  Of  particular  interest  is  the  behavior  of 
the  fundamental  lag  frequency  of  the  blade,  since  the  lag  motion  is  very  important  in  blade 
dynamic  study.  As  the  size  of  the  damage  (%)  increases,  the  irequency  of  the  first  lag  mode 
is  increased.  This  is  thought  to  be  caused  by  changes  in  the  centrifugal  force  due  to  blade 
mass/stiffness  effects.  The  torsional  mode  shap>e  of  the  20%  tip  damaged  blade  (Damage  2) 
is  quite  different  from  that  of  baseline  and  Damage  1  (see  Figure  12).  This  is  because  of  a 
large  coupling  effect  between  the  second  lag  and  the  first  torsional  modes. 


Most  of  the  comprehensive  aeroelastic  codes  developed  to  calculate  helicopter  flight 
dynamic  response  and  stability  employ  the  assumption  that  the  response  of  dll  the  blades 
is  identical  with  an  appropriate  phase  shift  for  each  blade.  Thus,  in  the  rotor  response  cal¬ 
culation,  a  set  of  coupled  flap- lag- torsion  equations  corresponding  to  a  single  blcide  is  used. 
In  the  study  of  ballistically  damaged  rotor  subsystems,  the  motion  of  each  blade  has  to  be 
represented  by  an  independent  set  of  equations,  and  their  response  is  calculated  individually. 
'Fhis  type  of  solution  approach  is  referred  as  a  “multi-blade  formulation.”  This  formulation 
also  drastically  increases  the  size  of  equations  to  be  solved.  For  example,  the  size  of  system 
equations  quadruples  for  a  four-bladed  rotor.  In  the  present  study,  a  multi-blade  formulation 
proposed  by  Chopra  [7]  is  used  to  calculate  the  rotor  hub  forces  and  moments.  Therefore, 
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the  motions  of  the  individual  blades  are  considered  in  the  following  discussion  of  results. 


Figures  13  through  15  present  the  time  histories  of  hub  loads  (nondimensionalized  by 
mQ^R^)  transferred  to  the  airframe.  Again,  for  comparison  purposes,  the  results  of  un¬ 
damaged  rotor  system  are  shown,  along  with  the  two  damage  cases  for  a  high  speed  flight 
condition  {ft  =  0.338):  Figures  13  and  14  for  the  in-plane  hub  loads  and  Figure  15  for  the 
out-of-plane  loads.  These  forces  were  calculated  using  a  multi-blade  formulation  considering 
the  motions  of  all  the  blades.  These  forces  cause  vehicle  vibration  and  are  very  important 
for  establishing  control  margins.  One  observes  the  3/rev  variation  of  forces  with  this  three- 
bladed  rotor  for  the  undamaged  (baseline)  case.  One  also  observes  a  significant  increase  in 
the  magnitude  of  all  three  force  components  as  the  amount  of  damage  increases  (Damage 
2  case).  It  is  also  interesting  to  note  that  the  vertical  force  variation  of  Damage  1  is  quite 
different  from  that  of  Damage  2,  as  shown  in  Figure  15.  This  points  out  the  importance  of 
nonlinear  aeroelastic  analysis  for  damaged  helicopter  rotor  subsystems:  One  may  not  simply 
apply  aeroelasiic  vulnerability  data  linearly  to  similar  damage  cases. 


V.  Conclusions 


Simulated  ballistic  damage  effects  on  helicopter  rotor  blade  response  and  rotor  hub 
loads  are  investigated  using  a  comprehensive  aeroelastic  analysis.  The  analysis  is  based  on 
finite  element  theory  in  space  and  time  coordinates.  Each  blade  is  treated  as  an  elastic 
beam  undergoing  combined  flap-lag-torsional  motions.  Multi-blade  formulation  is  used  to 
calculate  the  forces  and  moments  acting  on  the  rotor  hub.  Results  are  calculated  for  the 
baseline  undamaged  rotor  system  of  SA  349/2  Gazelle  helicopter  and  with  two  simulated 
blade  damage  conditions.  Ba-sed  on  the  analysis  results,  the  following  conclusions  are  drawn 
from  this  study: 


1.  Blade  damage  affects  all  the  modes  of  the  blade.  The  effects  are  more  distinct  in  the 
higher  modes. 

2.  As  the  amount  of  blade  damage  increases,  the  fundamental  frequency  of  the  lag  mode  is 
increased. 

3.  Blade  dissimilarity  induces  a  large  l/rev  variation  in  the  hub  loads. 

Recommendation  for  future  study 


From  this  study,  it  was  shown  that  a  20%  rotor  blade  tip  damage  would  induce  a  severe 
vibration  in  the  study  aircraft  primarily  because  of  large  rotor  blade  dissimilarity.  This  will 
most  likely  limit  the  helicopter’s  flight  performance  in  terms  of  the  vehicle’s  mobility.  This 
severe  vibration,  however,  will  also  affect  the  aircraft’s  dynamic  stability  and  control.  Future 
study  of  investigating  these  phenomena  would  be  of  great  interest. 
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Figure  3:  Correlation  of  SA349/2  Shaft  Power. 
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Figure  4:  Correlation  of  Flt^  Response  -  Coning  Angle  /3o. 
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Longitudinal  Disk  Tilt  (deg)  Lateral  Disk  Tilt  (deg) 


Figure  6:  Correlation  of  Flap  Response  -  Longitudinal  Disk  Tilt  0ic- 
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10  %  Damage  Case 


a)  Damage  1 


b)  Damage  2 


Figure  7:  Schematic  of  Blade  Damage  Cases. 
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Rotating  Blade  Mode  Shapes:  First  Lag 


Figure  10;  First  Lag  Mode  Shapea  of  Undamaged  and  Damaged  Blades. 
Rotating  Blade  Mode  Shapes:  Second  Lag 


Figure  11:  Second  Lag  Mode  Shapes  of  Undamaged  and  Damaged  Blades. 
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Longitudinal  Hub  Loads  Torsional  Deflections,  rad 


Rotating  Blade  Mode  Shapes:  First  Torsion 


Figure  12;  Torsional  Mode  Shapes  of  Undamaged  and  Damaged  Blades. 
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Figure  13:  Variations  of  Longitudinal  In-Plane  Hub  Loads  (/i  =  0.338). 
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Vertical  Hub  Loads  F-y  :  Lateral  Hub  Loads 


Figure  14:  Variations  of  Lateral  In-Plane  Hub  Loads  {fi  =  0.338). 


Figure  15:  Variations  of  Out-of-Plane  Hub  Loads  (ft  =  0.338). 
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Appendix  A 


UMARC  Input  Data  Set 


Flight  Condition  Inputs 
Advance  ratio,  /i. 

Rotor  speed  ratio,  Q/ilo  :  Oo  is  a  nominal  rotor  rpm. 

Density  ratio,  p/p,  :  p,  is  a  standard  sea-level  air  density. 

Flight  angle.  Off  (rad). 

Tip  Mach  number,  Mup. 

Main  Rotor  Overall  Configuration  Inputs 

Rotor  type  (articulated  or  hingeless  or  bearingless). 

Number  of  blades,  N^. 

Solidity,  <t 
Lock  number,  7. 

Rotor  thrust  ratio,  CtI<^- 

Hub  location  (x,y,z)  with  respect  to  aircraft  c.g. 

Main  Rotor  Dlade  Configuration  Inputs 

Linear  twist,  di^  (rad):  positive  for  leading  edge  up. 

Precone,  /3p  (rad):  negative  for  droop. 

Root  cutout,  Rroot- 

Number  of  spatial  beam  elements,  n. 

Blade  Element  Structural  Properties  Inputs  :  required  for  each  n  elements. 

Element  length,  U. 

Chord  length,  c. 

Element  mass,  m. 

Chord-wise  offset  of  tensile  axis  ahead  of  elastic  axis,  £„ . 
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Chord- wise  offset  of  aerodynamic  center  behind  elastic  axis,  ej. 

Chord-wise  offset  of  blade  center  of  mass  idiead  of  elastic  axis,  eg. 

Axial  stiffness,  EA. 

Section  stiffness  constants,  EBi,  EB^. 

Warping  rigidity,  ECi. 

Section  warping  constant,  ECj. 

Flap-wise  bending  stiffness,  Ely. 

Lag  (Chord-wise)  bending  stiffness,  Ely. 

Torsional  stiffness,  GJ. 

Radius  of  gyration  of  blade  cross  section,  ib^- 
Dlade  cross-sectional  mass  radius  gyration,  km- 
Blade  cross-sectional  mass  radius  of  gyration  in  flap  direction,  ibm,- 
Blade  crosvsectional  mass  radius  of  gyration  in  lag  direction,  km^. 

Blade  sweep  angle,  6,u,  (rad):  negative  for  swept  back. 

Blade  droop  angle,  (rad):  negative  for  droop  down. 

Tail  Rotor  Configuration  Inputs 

Tail  rotor  radius  ratio,  Rtr/R  :  /2  is  a  main  rotor  radius. 

Tail  rotor  gear  ratio,  :  0  is  a  main  rotor  rpm. 

Tail  rotor  blade  twist,  (rad):  positive  for  leading  edge  up. 

Number  of  blades,  Nt,. 

Solidity,  (Tj. 

Taul  rotor  hub  location  (x,y,z)  with  respect  to  aircraft  c.g. 

Liean  lift  curve  slope,  ci,,.. 

Mean  chord  length,  C(r- 

Horizontal  Tail  Configuration  Inputs 

Horizontal  tail  area  ratio,  Sht/A:  /I  is  a  main  rotor  disk  area  (x/2®). 
Longitudinal  distance  between  aircraft  c.g.  and  tail  aerodynamic  center,  zsi. 
Mean  lift  curve  slope,  ci^,. 
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Mean  chord  length,  Ch|. 


Fuselage  Configuration  Inputs 

Effective  fuselage  mass,  Mj . 

Fuselage  roll  inertia,  /^. 

Fuselage  pitch  inertia,  la- 
Fuselage  lift  coefficient,  Cij . 

Fuselage  drag  coefficient,  Co,- 
Fuselage  parasite  drag  area,  f /A. 
Fuselage  side  force  coefficient,  Cy^- 
Fuselage  pitch  moment  coefficient. 
Fuselage  roll  moment  coefficient,  Ctj. 
Fuselage  yaw  moment  coefficient,  C„j. 
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INTENTIONALLY  LEFT  BLANK. 
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Appendix  B 


Sample  Gaxelle  Input  Data 

title  :  gazelle  rotor  (  sa349/2  helicopter  ) 

1  us  aero  5  fuselage  dof  dynamic  inflow 

1  (stability  analysis) 

*_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 

rotor_properties 

!  rotor_type  t  articulated  flap_lag_delta3 

1  rotor_n£une  ;  gazelle 

!  no_blade8  *  3  solidity  *  0.064 

I  lock_no  *  5.5  ct/sigma  ■  0.065 

I  cg_below_hub  ■  0.25 

1  cg_hub_of fset_x  ■  0.017  cg_hub_offset_y  »  0.0 


b 1 ade_p roper t ie s 

!  twist (deg)  *  0. 

!  lag_spring_con8tant  ■  0.0663 
!  delta3_con8tant  ■  0.0294 

!  root_cut(%)  «  10.0 


precone (deg)  ■  0.0 
lag_daii^er_constant  ■  0.0151 
pitch_8pring_con8tant  «  0.0150 
no_space_element8  ■  8 


element_length 

« 

• 

0.100000 

0.101500 

chord 

« 

« 

0.0667 

0.0667 

j 

eiy 

« 

0.000781 

• 

0.000597 

eiz 

• 

« 

0.088477 

• 

0.039430 

gj 

• 

• 

0.000724 

• 

0.001077 

ea 

• 

50.0 

• 

• 

50.0 

offset_cg&ea 

• 

• 

0.023 

-0.029 

1 

of fset_taiea 

• 

• 

0.0  0 

1 

of fset_ac&ea 

• 

• 

-0.060 

!  0.057 


0.100000 

0.176000 

0.101500 

0.212400 

0.118100 

0.069500 

0.0667 

0.0667 

0.0667 

0.0667 

0.0667 

0.0667 

0.000494 

0.000620 

0.000608 

0.000012 

0.003862 

0.051004 

0.026000 

0.031780 

0.039430 

0.045540 

0.017200 

0.242910 

0.001800 

0.001184 

0.001077 

0.000964 

0.052152 

0.070430 

50.0 

50.0 

50.0 

50.0 

50.0 

50.0 

0.003  -0. 

023  -0.026 

0.056  0. 

0  0.0 

0.0  0.0 

0.0  0.0  0 

O 

o 

o 

0.060  0. 

057  0.057 

0.010  0. 

0  0.0 

27 


1  sq_kin2 
! 


0.000381 

0.000248 


0.000180 

0.000249 


0.000248 

0.000155 


0.000248 

0.000544 


1  sq_kinl 
! 


O.OOOOllC 

0.0000030 


0.0000040 

0.0000170 


0.0000030 

0.0000001 


0.0000030 

0.0001090 


1  pl_chord_length  ■  0.01 

!  pl_axial_length  ■  0.01 


pl_height  *  0.02 


* 


fuselage_properties 
!  no  hub  dof  *  5 


1  body_seq 


1  2  3  4  5 


I  ef fective_x_maa8  ■  86.1 

I  effective_y_ma8fl  ■  86.1 

!  effective_z_maa8  »  86.1 

1  fuselagejpitch_inertia  »  10.2709 

1  fuselage_roll_inertia  ■  1.7681 


* 


tail_propertie8 

gear_ratio_tail_rot 

m 

15.2 

solidity_tail_rot 

m 

0.4644 

rad_ratio_tail_rot 

m 

0.0662 

cg_tail_rot_of fset_x 

m 

1.1162 

t a i l_rot_above_cg 

m 

0.0933 

twist ( deg ) _tai l_rot 

m 

0.0 

cl_tail_rot 

■ 

5.7 

area_ratio_tail_hor 

. 

0.0 

cg_tail_hor_off8et_x 

■ 

1.0 

cl_tail_hor 

m 

4.5 

cO  tail  hor 

m 

0.0 

* 


airfoil_propertie8 
1  airfoil  «  myairfoil 

I  table_look_up  :  on  no_air foils  -  1 

!  air foil_8 tart ing_location8  ■  .0 
!  cl_table_naine8  test 
1  cd_tabl9_naines  test 
t  cm  table  names  test 
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€.00 


!  ref_lift_curve_8lope 
!  le  recov  fact  ■  0.97 


9 


f 1 ight_condit ion 

I  advance_ratio  -  0.338  den8ity_ratio  ■  1.0 

1  rotor_8peed_ratio  ■  1.0 

1  tip_mach  ■  0.65  flight_angle  ■  0.0 

* _ _ _ _ _ _ 

aerodynamic 8 

*  uniform_inf low 

*  dree8_linear_inf low 

!  r ever 8e_f low 

1  l_circ_us 

*  l_ua_drag 

*  l_usjpm 

*  l_impul_u8 

*  nl_u8_te_8ep 

1  dyn_stall 

*  sing_peak_pre8_wake 

!  sing_peak_free_wake 

*  dual_j5eak_pre8_wake 

*  dual_peak_free_wake 

1  free  wake  iter  no  «  9 


trim_analy8i8 

*  coupled_rb 

1  nl_8truct_aero 
1  fir8t_flap_freq  ■  1.016 


29 


*  coupled_trim_scheitMi 
1  coupled_triin_schem*  • 
t  delta_control8(%) 

I  conv_crit 
!  re8p_conv_crit 
1  triin_conv_erit 

!  no_max_iter 
!  qs_aero_dainp 
!  us_aero_dainp 

!  no_time_elemant8 
!  nodes_per_time_elem 

!  no_flap_mode8  ■  2 

1  f lap_mode_8eq  «  2  3 

t  no_lag_jnode8  *  2 

1  lag_mode_8eq  ■  1  6 

!  no__tor8ion_mode8  ■  1 
!  tor8ion_roode_8eq  ■  4 

1  no_axial_ntode8  -  0 
I  axial_mode_8eq  ■  10 


* 


stability_analysi8 

l  £ix_qs_const 
!  no_j)8i_location8  ■  24 
1  dynamic_inf low 


neifton 

marq_n«%^on 

5 

activ* 

0.005 

0.0001 

50 

40.00 

40.00 

8 

6 


1 


output_option8 


pr int_cont rol 
hub  load  control 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria.  VA  22304-614S 

1  Comnumder 

U.S.  Army  Materiel  Command 
ATTN:  AMCAM 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

1  Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-AD. 

Tech  Publishing 
2800  Powder  Mill  Rd. 

Adelphi,  MD  20783-1145 

1  Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-AD. 

Records  Management 
2800  Powder  Mill  Rd. 

Adelphi.  MD  20783-1145 

2  Commander 

U.S.  Army  Armament  Research, 
Devdopment,  and  Engineering  Center 
ATTN:  SMCAR-IMI-I 
Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Benet  Weapons  Laboratory 
U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCB-TL 
Watcrvliet.  NY  12189-4050 

1  Director 

U.S.  Army  Advanced  Systems  Research 
and  Analysis  Office  (ATCOM) 

ATTN:  AMSAT-R-NR,  M/S  219-1 
Ames  Research  Center 
Moffett  Field,  CA  94035-1000 


1  Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOO 
Redstone  Arsenal.  AL  35898-5010 

1  Commander 

U.S.  Army  Tank-Automotive  Command 
ATTN:  AMSTA-JSK  (Armor  Eng.  Br.) 
Warren.  MI  48397-5000 


1  Director 

U.S.  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSR 

White  Sands  Missile  Range,  NM  88(X)2-5502 

(cum-  •■ty)  1  Commandant 

U.S.  Army  Infaittry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 

Fort  Benning,  GA  31905-5t^ 

(UncUo.  only)  j  Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSH-WCB-O 
Fort  Benning,  GA  31905-5000 


1  WL/MNOl 
EglinAFB.FL  32542-5000 

Aberdeen  Proving  Ground 

2  Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

1  Cdr,  USATECOM 
ATTN:  AMSTE-TC 

1  Dir.  ERDEC 

ATTN:  SCBRD-RT 

1  Cdr,  CBDA 

ATTN:  AMSCB  ClI 

1  Dir.  USARL 

ATTN:  AMSRL-SL-I 

10  Dir.  USARL 

ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 
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No.  of 
Copies 

3 


2 


1 


1 


1 

2 


1 

1 


Organization 

Commander 

U.S.  Army  Aviation  Systems  Command 
ATTN:  AMSAV-ESC,  G.  Kovacs 
SFAE-AV-BH-T.  R.  Olson 
SFAE-AV-AAH-SA.  D.  Roby 
4300  Goodfellow  Blvd. 

St  Louis.  MO  63120-1798 

Commander 

AATD 

ATTN:  AMSAT-R-TF,  G.  Hufsietfcr 
SAVRT-TY-ASV,  H.  Reddick 
USAATCOM 

Fort  Eustis,  VA  23604-5577 

Commander 

ASC 

ATTN:  ASC/XRM,  G.  Bennett 
AFSC  Headquarters 

Wright-Patterson  AFB,  OH  45433-6503 

Commander 

WRDC 

ATTN:  WL/HVST.  M.  Lentz 
AFSC  Headquarters 

Wright-Patterson  AFB,  OH  45433-6553 
Commander 

ATTN:  ASC/XRYA.  K.  McArdle 
EglinAFB.FL  32542 


No.  of 

Copies  Ootanizmion 

1  Sikorsky  Aircraft 
Mail  Stop  ZIOIA 
ATTN:  G.  Burblis 
6900  Main  St 
Stratford.  CT  06601-1381 

Aberdeen  Proving  Ground 

1  Dir.  USAMSAA 

ATTN:  AMXSY-AD,  C.  Alston 

43  Dir.  USARL 

ATTN:  AMSRL-SL-B  (2  qjs) 

(Bldg.  328) 

AMSRL-SL-BA  (20  cps) 

(Bldg.  1065) 

AMSRL-SL-BA.  K.  Kim  (20  cps) 
(Bldg.  1065) 

AMSRL-SL-N.  M.  Miller 
(Bldg.  E3331) 


Commander 
Dahlgren  Diviskm 
Naval  Surface  Warfare  Center 
ATTN:  CodeG-13. 

T.  Wasmund 
D.  Dickinson 
Dahlgren,  VA  22448 

McDonnell  Douglas  Helicopter  Company 
ATTN:  William  Sims 
5000  E.  McDowell  Rd. 

Mesa,AZ  85205 

Boeing  Helicoptm  Division 
ATTN:  N.  Caravasos 
P.O.  Box  16858 
Philadelphia.  PA  19142 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-235 _ Date  of  Report  October  1993 _ 

2.  Date  Report  Received _ 

3.  Docs  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State.  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 

Organization 

OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 

(DO  NOT  STAPLE) 
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ATTN:  AMSRL-OP-CI-B  (Tech  Ub) 
Aberdeen  Proving  Ground,  MD  21005-5066 


